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The scattering effects of charm mesons by hadrons such as the pion, rho meson and 
nucleon are studied in an effective Lagrangian, and are found to be important in 
heavy ion collisions. At the CERN-SPS energies, hadronic rescatterings are shown 
to harden the charm meson mj< spectra, leading to a significant enhancement of 
the yield of intermediate-mass dimuons from charm meson decays. 



1 Introduction 

Heavy quark production in hadronic reactions is reasonably well described by 
perturbative QCD. However, in heavy ion collisions, final-state interactions 
may affect the spectra of produced heavy mesons. At the Relativistic Heavy 
Ion Collider (RHIC), a dense partonic system, often called the quark gluon 
plasma (QGP), is expected to be formed at the early stage. Since the QGP may 
induce a strong radiative energy loss of the produced heavy quarks, a change 
in the spectra of heavy meson observables could provide us information on 
the properties of the QGP. But interactions between heavy mesons and other 
hadrons may not be negligible and need to be studied. 

Such a hadronic modification of the charm spectra in heavy ion collisions 
has been recently suggested H as a possible explanation for the observed en- 
hancement of dimuons of intermediate masses in the NA50 experiments at 
the CERN-SPS. Assuming that charm mesons develop a transverse flow due 
to rescatterings with hadrons, leading to a harder charm meson mi spectra, 
dimuons from charm meson decays are also found to have a harder p± spec- 
trum. Based on the energy cuts for muons at the NA50 experiment, more 
dimuons would then be found to have an invariant mass above 1.5 GeV. An- 
other explanation based on dilepton .productions from secondary meson- meson 
interactions has also been proposed^]. Whether or not charm mesons acquire 
a transverse flow depends on how strongly charm mesons interact with other 
hadrons during their propagation through the matter. In this study, we shall 
evaluate the scattering cross sections of charm mesons with pion, rho, and nu- 
cleon, using an effective Lagrangian. The effects of hadronic scatterings on the 
charm meson transverse momentum spectra and dileptons from charm meson 
decays are then estimated for heavy ion collisions at CERN-SPS energies. 
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2 Charm meson interactions with hadrons 



We consider the scattering of charm mesons (D + , D~ , D°, D°, D* + , D*~ , 
D*°, and D*°) with pion, rho, and nucleon. If the SU(4) symmetry were exact, 
interactions involving pseudo-scalar and vector mesons could be described by 
the following Lagrangian: 



Cpp V = igTr (P^V^dpP) + h.c. 



(1) 



where P and V represent, respectively, the following 4x4 pseudo-scalar and 
vector meson matrices: 
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Expanding the Lagrangian in Eq.(|l]) in terms of the meson fields explicitly, 
we obtain the following Lagragians for meson-meson interactions: 



£ttDD* 
£ P DD 

r 



[D{d^) - (dpD)Tf] + h.c. , 



ig pDD [Df{dpD) - (dpD) TD]-p* , 



9 pirn P^ 



We also need the following Lagrangians for meson-baryon interactions: 
C-kNN = —ig^NNN^rN ■ n , 



£ p nn = 9 p nnN \^ff • pp + t^v^t ■ dpp u j N . 

Fig.[j] shows the Feynman diagrams considered in this study for charm 
meson interactions with the pion (diagrams 1 to 8), the rho meson (diagrams 
9 and 10), and the nucleon (diagrams 11 to 13)-. The differential cross sections 
of these processes can be found in the referenceB. Possible interferences among 
diagrams 3, 4 and 8 have not been included. 




Figure 1: The included diagrams for (a) Dir, (b) Dp, and (c) DN scatterings. 

For coupling constants, we take g p7r7r = 6.1, c/nD d* = 4.4, g P DD = 2.80, 
9ttNN = 13.5, g P NN — 3.25, and n p — 6.1. The SU(A) symmetry assumed in 
the Lagrangian in Eq.(Q) would give the following relations: 

9„KK*(3.3)=g„DD*(4A)=g pKK (3.0)=g p nrt2-8) = ^(3.4) = ^(3.0). 

The empirical values given in parentheses agree reasonably well with this pre- 
diction even though the SU(4) symmetry is badly broken. Form factors are 
introduced at the vertices to take into account the structure of hadrons. For t- 
channel vertices, monopole form factors of the form f(t) = (A 2 — m 2 )/(A 2 — t) 
are used, where A is a cut-off parameter a, and m a is the mass of the ex- 
changed meson. It should be noted that the cross sections for diagrams 2 and 
9 (D*tt <-> Dp) are singular because the intermediate mesons can be on-shell. 
In the present study, we simply add an imaginary part of 50 MeV to the mass 
of the intermediate pion as the regulator. 

The thermal averaged cross section, <av>, is shown in Fig.||(a) for initial 
particles with a thermal distribution at temperature T. Only the dominant 
scattering channels which have values above 1.1 mb are shown. 



3 Estimates of Rescattering Effects 

In this section, we estimate the effects of hadronic rescatterings on both the 
charm meson m± spectra and the invariant mass distribution of dileptons from 
charm meson decays in heavy ion collisions at SPS energies. 

We first determine the squared momentum transfer to a charm meson, p 2 ,, 
as the squared momentum of the final charm meson D2 in the rest frame of 
Di for a scattering process D\X\ — > D^X-i. In the charm meson local frame, 
we assume the time evolution of the hadron densities as p(r) oc 1/r. Then the 
total number of scatterings for a charm meson is given by 



N = I avpdr ~ avpoTo In 

J TO 
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Figure 2: Thermal averages (a) <o-t;>, and (b) < crvpg > for the dominant charm meson 
scattering processes as a function of temperature. The numbers in the notation of the curves 
correspond to the diagram numbers in Figure M. 



and the squared total momentum transfer from hadronic scatterings is 

f Ra < > 



< pi >=< Npl >-- 



= TT,p,N--- 



r In 



V T <P I 1> 



Thus, the relevant quantity is the thermal average < ovp^ > instead of the 
usual <<rv>. Fig.||(b) shows this thermal average for the dominant scattering 
channels which have values above 0.75 mb-GeV 2 . Summing up contributions 
from all scattering channels in Fig|l] (a), (b) and (c) separately, and simply 
dividing by 2 to account for the average over D and D* , we get 

< av p 2 >~ 1.1, 1.5 and 2.7 mb • GeV 2 at T = 150 MeV 

for 7r, p and N scatterings with charm mesons, respectively. 

For central Pb + Pb collisions at SPS energies, the initial total num- 
bers of particles are about 500^r), 220(p), 100(w), 80(7?), 180(iV), 60(A), and 
130(higher baryon resonances) U. We have not calculated the scattering cross 
sections between charm mesons and hadrons such as kaons, u), 77, A, and higher 
baryon resonances. For a conservative estimate of the effect, we only include 
7r, p and nucleon, and we obtain 

p Q T ~ 0.79(7r), 0.35(p), and 0.28(nucleon) fin" 2 , 
^< p 2 s >~ 0.61 GeV 2 =*> T s = 96 MeV . 

In the above, we have taken tq = 1 fm. The parameter Tg characterizes the 
scattering strength and is given by T$ ~< Pg > /(3tod) in the lowest-order 
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approximation 13. Based on Monte Carlo simulations, Ts has been related to 
the inverse slope T c g of the final charm meson m± spectrum, and this is shown 
in Fig. 6 of Ref.EJ. From that figure, we find that the charm meson T c g increases 
from 160 MeV to about 235 MeV, leading to a dimuon enhancement factor of 
about 2.1 for the NA50 acceptance. 

For heavy ion collisions at RHIC energies, in addition to hadronic rescat- 
terings of charm mesons, partonic rescattering effects on charm quarks also 
need to be included. Furthermore, radiative processes of charm quarks inside 
the QGP would further complicate the issue as they may cause energy loss and 
soften the charm meson m± spectra. More studies are thus needed before one 
can make predictions for RHIC. 



4 Discussions and Summary 

From our calculated scattering cross sections of charm mesons with hadrons 
such as pion, rho meson and nucleon, we have given an estimate of the rescat- 
tering effect at the SPS energies. We find that hadronic rescatterings in heavy 
ion collisions are likely to have a significant effect on charm meson spectra, 
and also the dilepton spectra from charm meson decays. 

The estimates given above are, however, based on a simple assumption on 
the time evolution of the hadronic system, which enables us to make a more 
analytical estimate for the rescattering effects. For a quantitative study of 
the rescattering effects on charm meson observables, studies with a partonic 
and hadronic cascade program are much needed as the time evolution and the 
chemical equilibration of the dense system are better simulated in such a model. 
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